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1. Introduction 

Confinement and spontaneous symmetry breaking as well as their respective counterparts - decon¬ 
finement and chiral restoration - are central issues of QCD. Their mutual relationship is among 
the main questions investigated in heavy ion collisions (HICs). A particularly challenging task 
is the establishment of the QCD phase diagram in a region of temperature T and baryo-chemical 
potential /./ which is accessible in HICs. While the beam energy scan at RHIC and forthcoming 
accelerator/detector installations, e.g. at FAIR/CBM and NICA/MPD, aim at finding signatures of 
the onset of deconfinement as a first-order phase transition in a critical point (CP), QCD based [1,2] 
and QCD rooted [3] calculations try to quantify its coordinates. Instead of facing the sign prob¬ 
lem of QCD at non-zero baryon density, one can resort to suitable models which display a CP 
at (7(p. jJ-cp) > (0,0) and a line of first-order phase transitions for /./ > /./( P (cf [4, 5] for recent 
surveys). By universality arguments, one can elucidate observables signaling the CP. 

In the previous search for deconfinement effects in HICs, electromagnetic probes (cf [6-8] for 
recent evaluations and further references) have been considered as promising since, by their pene¬ 
trating nature, they monitor the full space-time evolution of matter in the course of HICs. Given the 
proximity of deconfinement and chiral restoration at fi/T <C 1, the dilepton spectra have been con¬ 
sidered in [9,10] as candidates for messengers of chiral restoration. It appeared that at the pseudo- 
critical temperature of T pc ~ 154MeV the dilepton emission rates of confined matter (hadrons) and 
deconhned matter (quarks and gluons) become degenerate - quite natural since T pc quantifies the 
cross over location at pt/T <C 1. This degeneracy can be named quark-hadron duality and has been 
employed [11-13] to verify a “thermal radiation” component besides hard initial yield and electro¬ 
magnetic hadron final-state decays after the disassembly of the fireball in HICs. The approximate 
degeneracy of real-photon emission rates of confined and deconfined matter has been pointed out 
quite early [14] and can be used analogously to arrive at a consistent modeling of thermal radiation 
in HICs, either schematically [15,16] or with many refinements [17-19]. 

The recent analysis in [17,20], however, revealed that the calculated photon spectra fall short in 
comparison with data. The tension is caused, to some extent, by the consistency requirement of 
photon-V2 and photon-/; L systematics. According to [17,18], a solution is offered by the hypothesis 
of a “pseudo-critical enhancement” of the photon emission at T pc ± AT with AT = f/( 15 McV). It 
is therefore tempting to look for arguments which could support the hypothesis launched there and, 
more general, to extend the consideration over the QCD phase diagram and to check whether the 
CP can have a distinguished impact on the photon emission rate. 

2. Linear a model 

We utilize here the linear sigma model (LaM), which has been originally designed to model spon¬ 
taneous chiral symmetry breaking [21] and has been later coupled to the fermionic sector. It is a 
special quark-meson model with the Lagrangian = y + rM + m ++ ^s.b .> where 

the kinetic terms T£ r ) x „ and 2z fg M arc supplemented by the 0(4) symmetric meson self-interaction 
2z4r, which refers to an SU (T)i x SU(2)r symmetry which, in combination with the symmetry 
breaker T£ s .b. = Ha, dynamically generates mass terms; the fermion-meson coupling of Yukawa 
type is encoded in 2z f^M- The field content of Jz? is a two component Fermion (light quark doublet 
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\j/), an iso-vector (n) and an iso-scalar (a) meson. (The model can be extended to SU(3) and may 
include vector mesons [22] in the spirit of the gauged LaM [23] or can be supplemented by some 
gluon dynamics encoded in the Polyakov loop [24] or a glueball condensate [25].) When including 
linearized meson field fluctuations, as in [26-28], the resulting phase diagram has, in fact, a line of 
first-order phase transitions (cf solid curves in Fig. [I]), where the pressure as a function of T and /./ 
has a self-intersection defining T c {fj.). The fine structure of the region, where the CP is located [28], 
is thought to reflect the restriction to linearized fluctuations. (For a better account of fluctuations 
within a FRG approach, cf. [29].) 

The parameters in the Lagrangian are commonly fixed by requiring a coincidence of the dynam¬ 
ically generated meson quasi-particle masses in the vacuum (T = [i = 0) with the experimental 
values of pion, m™ c = 138MeV, and sigma, m™ c = 700MeV, and the constituent quark mass, 
adjusted to one third of the nucleon mass, = 312MeV. The remaining free parameter is deter¬ 
mined by (a)vac —fn- With these settings, the CP coordinates are {Tcp.pcp) ~ (74MeV,278MeV) 
(without the Fermion vacuum loop; for its impact cf. [30]). Since the sigma field in .5£\_o\\ can be 
considered as an effective field which may be tentatively identified with /o(600) leading to vacuum 
mass values of m™ c = 400 — 800MeV. It is instructive to study its impact on the phase diagram 
{cf. left panel of Fig. []]). It turns out that the phase contour curves are shifted approximately lin¬ 
early with m™ c . For an estimate we note A T pc ~ 0.1 2Am a on the temperature axis. The CP is 
shifted approximately parallel to the /.( axis by about 40 MeV per 100 MeV change of mf iC in this 
mass region. The impact of m f c is discussed extensively in [27]. We refrain here from a discussion 
of the phase structure under variations of all parameters at fixed r pc , which we define here as the 
curve, where the normalized heat capacity as a function of T at fixed /./ has a maximum - other 
definitions of the pseudo-critical temperature are conceivable. For further discussion of the phase 
diagram within effective chiral quark-meson models cf. [31-33]. 

Analogously to the mean field approximation [34] a crucial role is played by dynamically gener¬ 
ated mass terms, which may be considered as quasi-particle masses c {cf. Fig. | ) of thermal 
excitations of the \j/, n and a fields (for details see [35]). Prominent features are (i) degeneracy of 
m n and m a at temperatures larger than r pc (ju) or T c (fi), respectively, (ii) a rapid dropping of 
at T = T c (ji) and a mild dropping at T ~ r pc (ju) and (iii) a global minimum of m a at (Tcp.pcp) 
(for the rationale of the softening of the a-mode, cf, e.g., [36-38] and for the impact on di-photon 
spectra [39,40]). 

It is remarkable that the mass ordering changes over the T-p plane, see Fig. [I] (right panel). The 
delineation of 2m CJ > m n a vs. 2 m q < m n a coincides with T c {ji). There is a nar row valley, wherein 
m a < 2 m q < m n separating 2m q > m n a vs. 2 m q < m n a regions. This and the other delineation 
curves have no direct relation to T c (ji) or T pc ( ji ). 

3. Photon emission rates 

In a phenomenological approach and as a first exploratory step one can consider the above men¬ 
tioned quasi-particle masses as masses of the quanta entering the respective dispersion relations 
Ey x a = p 2 + ml, n a , which in turn show up in the thermal occupation functions (cf. appendix B 
in [29]). Due to the strong variations of these masses over the phase diagram {cf. Fig. |2j), one can 
expect also a non-trivial pattern of the real-photon emission rate (see figure 2 in [35]). 
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Figure 1: Left panel: Phase structure of the LcjM for various values of the sigma mass parameter m™ c = 
400, 500, 600, 700 and 800 MeV (bottom to top). The solid curves depict first-order phase transition lines 
and the dotted curves are an estimate for a pseudo-critical temperature T pc (p). Both curves meet at the CP 
(dots). Right panel: Relations of the effective masses m VJI ^ a over the phase plane (for m™ c = 700MeV, 
with phase sUucture as in left panel). 



Figure 2: Quark (blue solid curve), sigma (red dashed curve) and pion (green dotted curve) effective masses 
as a function of the quark chemical potential p for temperatures slightly below (left panel), at (middle panel) 
and slightly above (right panel) the critical temperature 7cp. The position of the critical chemical potential 
pep is depicted by the thin vertical line in the middle panel. 


Leaving a theoretically sound approach from Wightman functions to imaginary part of the retarded 
photon propagator to kinetic theory expressions from cutting rules for separate work, we employ 
the phenomenologically anticipated kinetic approach by estimating the rates via 

° Vjfc 3 =C l ^^^MP^f2(P2)a±f3(P3))\ M 12^y\ 2 S {4 \ Pl +P2-p3- k )A^ 

with C = (2(27T) 8 ) ', i.e. we account for lowest-order tree-level 2 —>• 2 processes 1 + 2 —>• 3 + y, 
where the photons are minimally coupled to J2 ?lctm according to [41], These are annihilations 
Yi + Yj —>• M + y and Compton processes Yi + Af —>• Yj + 7 and Vi + M —> Yj + 7- respectively, 
with i,j £ {u.d} and M € {tt°, k + .o) being the respective meson(s) fulfilling charge conservation 
depending on the choice of i,j. (Clearly, this list of reactions is by far not exhaustive when having 
in mind hadronic sources in general. For instance, in [42] the non-linear a model with vector and 
axial-vector mesons has been used to investigate real-photon emission rates from hadron sources. 
In [43], the channel n + 7i~ —» o /p —>• K ] n y has been identified as important for soft photons. This 
and many more channels should be considered when attempting a comprehensive rate estimate.) 
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The matrix elements for Compton and annihilation processes are related by crossing symmetries. 
The evaluation of ( |3.l| ) by means of the S (4] and exploiting the cylinder symmetry in the center of 
mass frame of the reaction leaves four integrals [44,45] to be done numerically. 

4. Photon spectra 

An instructive form of ( |3.1| ) can be derived by resorting to the Boltzmann approximation: 

( 0 ^ « Ce^' T [~ds^^e-W 7 (4.1) 

CIK J sq S JTl^ 

withC = T/(32(0(2k) 6 ), c(s ) = //+ dt\M[ 2 -> 3 y\ 2 , t± = m\ — s -^-[(s+m\—m^) + \J X(s,m 2 ,m%)J, 
X(x,y,z) = x 2 +y 2 + z 2 — 2xy — 2xz-2yz, F(s ) = (s — m 2 ) / (A(o) + (s( 0 ) / (s — m 2 ) and y/so = 
max{mi + n? 2 ./?? 3 }. The vaiiable e discriminates between three reaction types under consider¬ 
ation: annihilation (e = 0), Compton scattering at quarks (e = +1) and Compton scattering at 
antiquarks (e = —1). For disentangling effects of the phase space distributions /i_ 3 , in particu¬ 
lar the role of the involved quasi-particle masses, and matrix elements A/ 12 —> 3 y encoded in a ( 5 ) 
let us analyze the function F. Two cases are to be distinguished: (i) If the position of the min¬ 
imum of F w.r.t. s, .V| = m 2 + 2 (Om 3 , lies within the range of integration, i.e.so < si, a Taylor 
expansion F(s) = F(s 1 ) + 1/2 • F"(s 1 ) • (s — sq) 2 + ff((s — jq) 3 ) is in order and the exponential 
in ( [4. lj ) can be replaced by a Gaussian with width \J 4m 3 ft ) 2 T and height exp{—(ta + m 3 )/T}, 
making the exponential thermal suppression apparent, (ii) If si < so, a Taylor expansion F(s) = 
F(so) + F'(sq) ■ (s — so) + ff((s — so) 2 ) yields a s independent suppression factor exp{— F(sq)/T} 
that develops a maximum w.r.t. (0 at (O p = (sq — m 2 ) /(2^/Tq), meaning the spectrum displays a peak 
at (O p > 0. Whether it appears is thus controlled by the relative size of so and ,sq. If m 1 +mj < m 3 , 
then sq = m 2 < si (cf case (i)), and the photon rate has its maximum at (0 = 0 and declines ex¬ 
ponentially with power law corrections originating from a(s) and other omitted factors in (p~l|). 
On the other hand, if m 1 +m 2 > m 3 , then si < so (cf. case (ii)), and the rate has a maximum 
oc e -( m i + m 2 )/T at Furthermore, even for ni\ +m 2 > m 3 , there will be a certain photon frequency 
(O s > (O p for which s () = .sq, implying that for (0 > (O s case (i) is applicable and the rate behaves as 
oc exp{— (o/T} at large ( 0 . Thus the mass ordering determines some basic features of the spectra, 
which is the reason why we have shown in the right panel of Fig. [I] the corresponding regions. 

In Fig. Rj. the photon spectra originating from four selected processes are plotted for three different 
positions in the phase diagram, corresponding to the chirally broken phase (a), the chirally restored 
phase (c) and the proximity of the CP (b). Comparing the spectra based on ([Q|), one notices that 
in fact some of the processes develop a maximum at ~ a> p and vanish at 0 ) 0. while others seem 

to diverge in the IR limit as anticipated above by generalizing a previous consideration in [45]. 
With the above approximations one can understand why at different positions in the phase diagram 
the IR behavior of the photon rates is so different. It is the interplay of the masses of the involved 
modes and to a less extent the details of the interaction process. Since for the Compton process 
the sum of the incoming masses m v + m nx: is larger than the mass of the outgoing quark m ¥ , the 
corresponding contribution to the photon rate always has a maximum at finite photon energies. For 
the annihilation process there is a sensible dependence on the position in the phase diagram: In the 
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Figure 3: Photon spectra co J ^ for the LetVI evaluated at Li = OMeV (a), p = 270MeV (b) and p = 
360MeV (c) for T = 80MeV; panel (d) is a zoom into (b). Depicted are annihilations into photon and sigma 
y/ + yr —» 7 + a (solid blue curves) and into photon and pions y/ + y/ —> y + n (red short dashed curves) as 
well as the corresponding Compton processes y/ + <y —► y/ + y (green long dashed curves) and y/ + Tt —> y/ + y 
(cyan dotted curves). 


chirally broken phase, the sigma mass can be larger than twice the quark mass (cf the diagonally 
hatched region in the right panel of Fig. |T[). Therefore, there is no maximum of the rate at non-zero 
photon energy. Instead, the exponential factor approaches some finite value at co = 0 and other - 
in the former case subleading - effects get dominant. The most prominent effect stems from IR 
divergencies of the matrix elements which let the photon rate diverge as co 2 at co —»■ 0 and thus 
needs to be regularized. In a band around the phase transition line and the pseudocritical region 
(horizontally hatched region in the right panel of Fig. [[]) the meson masses are less than twice 
the quark mass ( e.g . m v ~ m a ~ m n ~ 200MeV at (T.p) ~ (7cp,/fcp))- Therefore, all processes 
lead separately to spectra which vanish at co ^ 0 and have a maximum of similar height and an 
exponential tail at large co. From such a consideration one can already conclude that a divergency 
will not appear in a region of the phase diagram where the sigma mass is smaller than twice the 
quark mass. In the high temperature phase the quarks get light and the mesons heavy ( cf. Fig |2|). 
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Figure 4: cr as a function of ,v for the re¬ 
actions exhibited in Fig. (same linestyle) at 
the same positions in the phase diagram (p = 
0,270 and 360 MeV from (a) to (c) at T = 
80MeV). The arrows depict the maxima of 
the integrand in ( |4.1[ ) for co = 50 MeV, where 
the contribution to the respective emissivities 
are largest. 


Therefore, the annihilation rates have no maxima in the unhatched region of Fig. [I] (right panel) 
and hence are not exponentially suppressed at co —> 0 leading - together with the IR divergency of 
the matrix element - to a large photon rate at small co. Thus the pattern of the rate behavior for 
the subprocesses shown in Fig. [| can be explained quite naturally. Comparing parameter sets with 
different vacuum sigma masses, the above considerations imply a strong change of the low energy 
photon rates over the phase diagram if we contrast the emissivity of the LaM with low vacuum 
sigma mass (m™ c 1 $ 2 m ™ c ) to the emissivity at higher vacuum sigma mass (m™ c > 2 m™ c ). 

To highlight the role played by the matrix elements we also show in Fig. j| the .v dependence of 
a. Again there are two types of functions: If m\ +mi < m 3 , a(s) diverges in the limit s —>• so like 
~ (.v — .vq) ', and conversely, if mi + m 2 > m 3 , then a(s) behaves like y/s — sq. The reason for the 
divergence in the former case are IR divergencies of the matrix elements, because if v > vo = m? 
the photon energy in the center of mass frame is small and thus the matrix elements are enhanced 
because the numerators in the propagators get small. 
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Figure 5: Contour plots of photon rates (B ^J; in units of MeV 2 for (0 = 50MeV over the phase diagram 
for the reactions (a): y/ + n —> y+ yr, (b): y/ + iff —> 7 + k, (c): y/ + a —> 7 + yr, (d): yf+yf -> 7 + a. The 
phase structure is as in Fig. [I] (right panel). 

5. Variation of photon rates over the phase diagram 


The photon rates at fixed photon frequency depend sensitively on the position in the phase diagram 
as shown in [35] for co = lOMeV. A rate maximum in the critical region for the a-involving 
Compton process was found. We focus here on higher values of 0 ). From Fig. 3(d)| one notices 
that up to w ~ 50MeV the a-involving Compton process (for which the available phase space 
around the critical point is enlarged, because of the small sigma mass) is the dominant channel 
and hence one can hope that a signal characteristic for the critical point can be obtained. For this 
frequency, the rates are depicted in Fig. ||. Several features are apparent. In the high temperature 
phase the Compton rates (a) and (c) are suppressed by many orders of magnitude relative to the 
corresponding annihilation processes (b) and (d). While the ^-involving Compton rate is largest 
in the crossover region, the a-involving Compton rate exhibits a global maximum in the critical 
region. In the low-temperature phase, the annihilation into a a meson and a photon is the strongest 
contributing process and shows the remnant of the structure seen in figure 2 in [35]. All of these 
observations can be explained with the above reasoning. The suppression of the Compton processes 
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is explained by the different approximation schemes necessary for the two types of processes. Since 
mi +n ?2 is for the Compton process always larger than m 3 = mi, and co = 50MeV is far below the 
peak frequency in the high temperature phase the exponential suppression factor exp{— F(sq)/T} 
is many orders of magnitude smaller than the one for the annihilation case exp{—(m 3 + to) /T }. 
For 0 ) = f/( \ — 2GeV) our results show no special features at T pc for small /./. 

6. Summary 

Employing the linear sigma model (LaM) we investigate whether the soft-photon emission rates 
can reflect the conjectured phase structure of QCD. The LaM is chosen as a simple approach which 
exhibits a critical point (CP) at non-zero chemical potential. Relying on the LaM field content, 
which is very schematic and mirrors only in a limited manner the proper degrees of freedom of 
QCD, we find, however, that the changes of the quasi-particle excitations masses within the phase 
diagram give rise to significant changes of the emission rates in selected channels. In particular, 
the spectral shapes depend strongly on the effective masses of excitation modes involved. This 
lets us hope that more advanced considerations can dig out peculiarities of the total emission rate, 
e.g. caused by the softening of the a-type mode near the CP. 

The work is supported by BMBF grant 05P12CRGH1. 
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